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Abstract 

We present preliminary measurements of B meson decays to B + — ► luit + , B + — > u;JC + , and B° — > 
wif . The data were recorded with the BABAR detector and correspond to 88.9 x 10 6 BB pairs 
produced in e + e~ annihilation at the T(4S) resonance. We find statistically significant signals for 
all three channels: B(B + -► ujtt + ) = (5.4±1.0±0.5) x 10~ 6 , B(B + -> uK + ) = (5.0±1.0±0.4) x 10~ 6 , 
and B(B° — > loK°) = (5.3]ti'2 ±0.5) x 10 -6 . We also measure time-integrated charge asymmetries 
Ach(B+ -» um + ) = 0.04 ± 0.17 ± 0.01 and A ch {B + -> uif+) = -0.05 ± 0.16 ± 0.01. 
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1 Introduction 



We report the results of searches for B decays to the charmless final states u>ir + , ujK + , and ujK°. 
We reconstruct the lo mesons via the dominant decay mode u — > 7r + 7r~7r°, and K° via K° — > 1^ — > 
7r + 7r~. For the charged modes we also measure the direct CP- violating time-integrated charge 
asymmetry A ch = {V" - r+)/(r _ + where T ± = T{B ± -► w/i*). 

Table 1 summarizes the current knowledge of these decays, coming from measurements by 
CLEO [1], BABAR [2, 3] and Belle [4]. CLEO and BABAR find significant signals for the B + -► wtt + 
channel. The B + — > u;i^ + decay has an interesting history. It was originally seen by CLEO but 
a re-analysis with the full CLEO dataset could not confirm the earlier large branching fraction 
measurement. BABAR confirmed the smaller result of the CLEO re-analysis (see Table 1) but Belle 
has now published a significant observation with a relatively large branching fraction. A more 
precise measurement would help settle the situation. BABAR reported a significant signal for the 
B° — ► ujK° decay at conferences in summer 2002. 

The early indications from CLEO and now Belle of a branching fraction for B + — > uK + of 
> 10 x 10~ 6 were hard to accommodate theoretically. The B + — > ujK + decay is an interesting 
penguin-dominated decay with cancellations between primary Wilson coefficients. The B + — > luit + 
decay is expected to be dominated by the external and color-supressed tree diagrams. Typical 
calculations found branching fractions for both oj decays of a few times 10" 6 [5, 6, 7], although 
parameter tuning was able to accommodate > 10 x 10 -6 with difficulty. The decay B° — > uiK° is 
expected to have a comparable branching fraction to that of B + — > uK + . However the presence 
of a (suppressed) tree diagram in the charged decay would reduce the rate than for the neutral 
decay if there is substantial destructive tree-penguin interference for B + — ► tuK + . The theoretical 
situation has changed little recently. There are few calculations for these modes from the QCD 
factorization or perturbative QCD groups. One recent calculation of pseudoscalar-vector modes 
with QCD factorization [8] finds predicted branching fractions of < 5 x 10~ 6 for all three decays. 

2 Detector and Data 

The results presented in this paper are based on data collected in 1999-2002 with the BABAR 
detector [9] at the PEP-II asymmetric-energy e + e~ collider [10] located at the Stanford Linear 
Accelerator Center. An integrated luminosity of 81.9 fb _1 , corresponding to 88.9 million BB pairs, 
was recorded at the T(4S) resonance ( "on-resonance" , center-of-mass energy y/s = 10.58 GeV). 
An additional 9.6 fb _1 were taken about 40 MeV below this energy ("off-resonance") for the study 
of continuum backgrounds in which a light or charm quark pair is produced instead of an T(45). 

The asymmetric beam configuration in the laboratory frame provides a boost of ((3^) = 0.56 
to the T(4S). Charged particles are detected and their momenta measured by a combination of 
a silicon vertex tracker (SVT), consisting of five layers of double-sided detectors, and a 40-layer 
central drift chamber, both operating in the 1.5 T magnetic field of a solenoid. Photons and 
electrons are detected by a CsI(Tl) electromagnetic calorimeter (EMC). 

Charged particle identification (PID) is provided primarily by an internally reflecting ring- 
imaging Cherenkov detector (DIRC) covering the central region, though the average energy loss 
(dE/dx) in the tracking devices is also used for the pion daughters from the u> decay. From the 
Cherenkov angle, a K-ir separation greater than 4 standard deviations (a) is achieved for tracks 
with momenta below 3 GeV/c, decreasing to 2.5 a at the highest momenta in the final states 
considered here [11]. 
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Table 1: Summary of branching fraction results for B decays to u> mesons from CLEO [1], previous 
BABAR measurements [2, 3], Belle [4], and the present analysis. The results for all fits are given 
as well as a 90% confidence level upper limit if the measured yield is not judged to be significant. 
The signal yields and efficiencies (e) are also given. 



Expt. 


# BB (xlO 6 ) 


Fit £(xl0- 6 ) 


UL ,B(xl0- 6 ) 


Signif. (a) 


Signal yield 


e (%) 


— ► LOK ^ 
















10 


O.Zl 1 g ± 0.O 




O 1 

z.l 


7 O+6-0 

7.9_ 47 


26 


dAdMx 


zo 


-i ^+1.3 _i_ n q 

1.4_^ g ± U.o 


Q Q 
O.O 


1.0 


a 4+5.6 


1 o 


Belle 


32 


9.2+23 ± 1-0 




6.0 


18 9 +5 A i 


6.0 


This result 


89 


5.0 ± 1.0 ± 0.4 




8.9 


87 ± 15 


18 


B+ OJ7T+ 














CLEO 


10 


11.3iy ±0.8 




6.2 


oq c+8.2 
Zo.O 7 3 


26 


BABAR 


23 


6.6^ ±0.7 




5.1 


28+jf 


19 


Belle 


32 


4.2+^ ± 0.5 


8.1 


3.3 


10-4111 


7.7 


This result 


89 


5.4 ± 1.0 ±0.5 




8.4 


101 ± 18 


19 


B° loK° 














CLEO 


10 


10.01^ ± 1.4 


21 


3.9 


7.0+1;° 


7.4 


BABAR 


62 


5.9±i;| ± 0.9 




6.6 


27 + f 


7.4 


This result 


89 


5.3l};^ ±0.5 




7.5 


33+^ 


7.0 



3 Event Selection 

Monte Carlo (MC) simulations [12] of the target decay modes and of continuum and BB back- 
grounds are used to establish the event selection criteria. The selection is designed to achieve high 
efficiency and retain sidebands sufficient to characterize the background for subsequent fitting. 

From photons with energy E 1 > 50 MeV, we select ir° candidates by requiring the invariant 
mass to satisfy 120 < m 77 < 150 MeV. Candidate to mesons must satisfy 735 < < 825 MeV. 

For Kg candidates we require 488 < < 508 MeV, the three-dimensional flight distance from 
the event primary vertex > 2 mm, and the two-dimensional angle between flight and momentum 
vectors < 40 mrad. 

We make several particle identification requirements to ensure the identity of the signal pions 
and kaons. Daughter tracks of uo candidates must have DIRC, dE/dx, and EMC responses con- 
sistent with pions. For the B + — > uK + decay, the prompt charged track must have an associated 
DIRC Cherenkov angle between —5a and +2 a of the value expected for a kaon. For B + — > u;7r + , 
the DIRC Cherenkov angle must be between —2a and ±5<r of the value expected for a pion. 

A B meson candidate is characterized kinematically by the energy-substituted mass rriES = 
yj (^s + p • pb) 2 /Eq — ~p B and energy difference AE = E* B — \\fs, where the subscripts and B 
refer to the initial T(45) and the B candidate, respectively, and the asterisk denotes the T(AS) 
frame. We require | AE\ < 0.2 GeV and 5.2 < m^s < 5.29 GeV. The resolutions on these quantities 
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are about 30 MeV and 3.0 MeV, respectively. 

3.1 Tau, QED, and continuum background 

To discriminate against tau-pair and two-photon background, we require that the event contain at 
least five (four) charged tracks for neutral (charged) B pairs. To reject continuum background, we 
define an angle 9t between the thrust axes of the B candidate and of the rest of the tracks and 
neutral clusters in the event, calculated in the center-of-mass frame. The distribution of | cos#t| 
is sharply peaked near 1.0 for combinations drawn from jet-like qq pairs, and nearly uniform for 
the isotropic B meson decays; we require | cos#t| < 0.8 for the charged modes and | cos#t| < 0.8 
for the lower-background B° — ► ojK° decay. A second B candidate satisfying the selection criteria 
is found in about 10-20% of the events. In this case the "best" combination is chosen as the one 
closest to the nominal oj mass. 

3.2 BB background 

We use MC simulations of B^B and B + B~ pair production and decay to look for possible BB 
backgrounds. From these studies we find that BB background is small for all three decay channels. 
The uncertainty is limited by the statistical errors of these studies and is included in the systematic 
errors (see §6). 

4 Maximum Likelihood Fit 

We use an unbinned, multivariate maximum-likelihood fit to extract signal yields for our modes. 
With the cuts in §3, candidates are selected to match the kinematic structure of the desired decay 
chain. 

4.1 Likelihood Function 

We incorporate several uncorrelated variables for the kinematics of the B decay chain and a Fisher 
discriminant T for the B production and energy flow. Thus the input observables are AE, mgs, 
miTTTn, T~t = | cos8h\, and a Fisher discriminant T. The helicity angle 6h is defined as the angle, 
measured in the oj rest frame, between the normal to the oj decay plane and the flight direction of 
the oj. The Fisher discriminant [13] combines four variables: the angles with respect to the beam 
axis, in the T(4S*) frame, of the B momentum and B thrust axis, and the zeroth and second angular 
moments Lq,2 of the energy flow about the B thrust axis. The moments are defined by 



where 9i is the angle with respect to the B thrust axis of track or neutral cluster i, pi is its 
momentum, and the sum excludes the B candidate. 

Since we measure the correlations among the observables in the data to be small, we take the 
probability density function (PDF) for each event to be a product of the PDFs for the separate 
observables. We define two hypotheses j, where j represents signal or continuum background. The 
product PDF (to be evaluated with the observable set for event i) is then given by 




(1) 



V) = V 3 (m ES ) • Vj(AE) • Vj{F) ■ V 3 (m^) ■ Vj(H) . 



(2) 
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The likelihood function for each decay chain is 



* 3 

where Yj is the yield of events of hypothesis j found by maximizing C and N is the number 
of events in the sample. The first factor takes into account the Poisson fluctuations in the total 
number of events. 



4.2 Fit Parameters 

The determination of PDF parameters for the likelihood fit is accomplished with use of Monte 
Carlo (MC) simulation for the signal and on-peak data from a tiies-AE sideband for the continuum 
background. Several of the principle background parameters are allowed to float in the final fit. 

Peaking distributions (signal masses, AE, and J 7 ) are parameterized with Gaussian functions, 
with a second or third Gaussian or asymmetric width as required for good fits to these samples. 
Because these are pseudoscalar-vector decays of the B, the signal helicity- angle distribution is 
proportional to 7i 2 . Slowly varying distributions (mass, energy or helicity-angle distributions for 
combinatoric background) have low order polynomial shapes, with the peaking and combinatoric 
uj mass components each having their own H shape. The combinatoric background in tues is 
described by a phase-space- motivated empirical function [14]. The background Fisher PDF contains 
a second Gaussian component wide enough to ensure the PDF is not excessively small anywhere 
in its domain. 

Control samples of B decays to charm final states of similar topology are used to test the 
quality of the MC simulation for variables describing B decay kinematics. Where the control data 
samples reveal differences from MC in mass or energy resolution, we scale the resolution used in 
the likelihood fits. 

The efficiency bias of the likelihood fit is determined from simulated data samples. The qq 
background in these samples is generated from the continuum PDF shapes. A small number of 
signal events from the detailed MC simulation are added to create a sample with the same size as 
the data. The bias in the fit efficiency is determined from the mean of the fit yield for 500 such 
simulated samples. 



5 Fit Results 

By generating and fitting samples with both signal and continuum background events generated 
from the PDFs, we verify that our fitting procedure is functioning properly. We find that the 
minimum — In C value in the on-resonance sample lies well within the — In C distribution from 
these simulated samples. 

In Table 2 we show the results of the fits for off-peak and on-peak data. Shown for each decay 
mode are the number of events that were fit, the signal yield, the fully corrected efficiency and 
product branching fraction (for the u> and K® decays), the measured branching fraction, and the 
statistical significance of the result. We also give the branching fraction after correction for a 
crossfeed between the two charged channels. The K/ir misidentification rate of 9 ± 2% is found 
in studies with kaon and pion samples tagged kinematically from the decay D* + — > ir + D° with 
D° — > K~ir + . The statistical error on the number of events is taken as the change in the central 
value when the quantity — 2ln£ changes by one unit. The statistical significance is taken as the 
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Table 2: Fit values for u>K + , cj7t + , and ujK°. The corrected B for the charged modes is the 
branching fraction after correcting for crossfeed from one charged mode into the other. 



Fit quantity 




UJTi + 


toK° 


rit sample size 








On-resonance 


16729 


30563 


9563 


Off-TPsnnanrp 

vy 11 -1 V_ OVJIJ-CIlIJ-V 


1900 


3490 


972 


Signal yield 






33±i 


On-res data 


87 ± 15 


101 ± 18 


Off-res data 


n n+ 2 - 7 
u - u -o.o 


n n+ 3 - 6 
u - u -o.o 


n n+ - 9 


Combinations /event 


1.13 


1.15 


1.15 


Selection e (%) 


20.0 


21.9 


23.1 


(%) 


88.8 


88.8 


30.5 


Corr. ex]JBi (%) 


17.8 


19.4 


7.0 


Stat, significance (a) 


8.9 


8.4 


7.5 


£(xKT 6 ) 


5.5 ± 1.0 ±0.4 


5.9 ± 1.0 ±0.5 




Corrected ^(xlO -6 ) 


5.0 ± 1.0 ±0.4 


5.4 ± 1.0 ±0.5 




Signal A ch 
Background A c h 


-0.05 ±0.16 ±0.01 
-0.005 ± 0.008 


0.04 ±0.17 ±0.01 
-0.012 ±0.006 





square root of the difference between the value of — 21n£ for zero signal and the value at its 
minimum. We also give the charge asymmetry A c h for the charged modes for both signal and qq 
background. 

In Fig. 1 we show projections of uies and AE made by selecting events with signal likelihood 
(computed without the variable shown in the figure) exceeding a mode-dependent threshold that 
optimizes the expected sensitivity. 

6 Systematic Uncertainties 

Most of the systematic errors on yields that arise from uncertainties in the values of the PDF 
parameters have already been incorporated into the overall statistical error, because their back- 
ground parameters are free in the fit. We determine the sensitivity to parameters of the signal PDF 
components by varying these within their uncertainties. The results are shown in the first row of 
Table 3. This is the only additive systematic error. 

The uncertainty in our knowledge of the efficiency is found from auxiliary studies to be 0.8AT t %, 
2.5iV 7 %, and 3% for a K® decay, where N t and iV 7 are the number of signal tracks and photons, 
respectively. We estimate the uncertainty in the number of produced BB pairs to be 1.1%. The 
systematic bias from the fitter itself (1-2%) is estimated from fits of simulated samples with varying 
background populations. Published world averages [15] provide the B daughter branching fraction 
uncertainties. We account for systematic effects in cos #t (1%) and in the PID requirements (0.5%). 
Values for each of these contributions are given in Table 3. 
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Figure 1: The B candidate uies and AE for B + — * ujtt + (a, b), B + — ► loK + (c, d), and B° — ► u>K° 
(e, f). Points with errors represent data passing a cut on a likelihood ratio calculated without the 
quantity that is shown in the plots. The solid curves show the projected fit functions and dashed 
curves the background functions. 

A study of the charge asymmetry as a function of momentum for all tracks in hadronic events 
bounds the tracking efficiency component of a charge-asymmetry bias to be below 1%. _D*-tagged 
D — ► Kit and B samples provide additional crosschecks that the bias is negligible. We assign a 
systematic uncertainty for A c h of 1.1% based on the tracking study and a small PID contribution 
determined from the D* studies. 

7 Conclusion 

We report preliminary measurements of branching fractions for the decays B + — > uir + , B + — > ujK + , 
and B° — > uK°. We find statistically significant signals for all three decays and measure the 
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Table 3: Estimates of systematic errors (in percent) for the uiK + , o;7r + and ujK° branching fractions. 



Quantity 




W7T + 




"Pit T^olr! 

r li yield 


9 7 
Z. I 


0. 1 


3 9 
O.Z 


Fit efficiency /bias 


9 1 

z. / 


o.y 


9 1 
Z.l 


Track multiplicity 


1 n 
l.U 


1 n 
l.U 


1 n 
l.U 


Tracking eff/qual 


z.4 


z.4 


o. f 


7T / 7 eff 


5.0 


5.0 


5.0 


Kg efficiency 






2.9 


IN umber BB 


1.1 


1.1 


1.1 


Branching fractions 


1.0 


1.0 


1.0 


MC statistics 


1.0 


1.0 


1.0 


cos #T 


1.0 


1.0 


1.0 


PID 


1.4 


1.4 


1.0 


BB Background 


1.1 


1.0 


3.0 


Total 


7.3 


8.0 


8.8 



following branching fractions: 

(5.4 ±1.0 ±0.5) x 10~ 6 , 
(5.0 ±1.0 ±0.4) x 10~ 6 , 
(5.3^ ±0.5) x 10~ 6 . 

These results supersede the previous BABAR measurements [2, 3]. The result for B(B + — ► ujK + ) 
is larger than the previous measurement with one quarter the luminosity. From studies of the old 
and new measurements, we conclude that the difference is a statistical fluctuation. 

These results are much more precise than previous results and in good agreement with the- 
oretical expectations. A branching fraction for B + — ► uK + near 10 x 10~ 6 , as reported first by 
CLEO and more recently by Belle, is now definitively ruled out. In addition, we measure the charge 
asymmetries A ch (B+ -► ujit + ) = 0.04 ± 0.17 ± 0.01 and A ch (B + -► ujK + ) = -0.05 ± 0.16 ± 0.01. 
These values are consistent with being small as generally expected for these decays [16] though 
both experimental and theoretical uncertainties are still large. 
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